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Abstract — In this paper, transmission over time-selective, flat 
fading relay channels is studied. It is assumed that channel fading 
coefflcients are not known a priori. Transmission takes place in 
two phases: network training phase and data transmission phase. 
In the training phase, pilot symbols are sent and the receivers em- 
ploy single-pilot MMSE estimation or noncausal Wiener filter to 
learn the channel. Amplify-and-Forward (AF) and Decode-and- 
Forward (DF) techniques are considered in the data transmission 
phase and achievable rate expressions are obtained. The training 
period, and data and training power allocations are jointly 
optimized by using the achievable rate expressions. Numerical 
results are obtained considering Gauss-Markov and lowpass 
fading models. Achievable rates are computed and energy-per-bit 
requirements are investigated. The optimal power distributions 
among pilot and data symbols are provided. 

I. Introduction 

In wireless communications, channel conditions vary ran- 
domly over time due to mobility and changing environment. 
If the channel conditions are not known a priori, practical 
wireless systems generally employ pilot symbols to learn 
the channel. In one of the early studies in this area. Cavers 
provided an analytical approach to the design of pilot-assisted 
transmissions in [1] and [2]. Considering adaptive coding of 
data symbols without feedback to the transmitter, Abou-Faycal 
et al. [4] studied the data rates achieved with pilot-symbol- 
assisted modulation (PSAM) over Gauss-Markov fading chan- 
nels. The authors in [6] also studied the PSAM over Gauss- 
Markov channels and analyzed the power allocated to the data 
symbols when the pilot symbol has fixed power. They showed 
that the power of a data symbol decreases as its distance to 
the pilot symbol increases. 

Recently, cooperative wireless communications has attracted 
much interest. Cooperative relay transmission techniques have 
been studied in [8] and [9] where Amplify-and-Forward (AF) 
and Decode-and-Forward (DF) models are considered. How- 
ever, most of the studies have assumed that the channel condi- 
tions are perfectly known at the receiver and/or transmitter. In 
one of the recent studies, Wang et al. [7] considered wireless 
sensory relay networks where the conditions of the channels 
are learnt imperfectly only by the relay nodes. 

In this paper, we study the training-based transmission and 
reception schemes over a priori unknown, Rayleigh fading 
relay channels in which the fading is modeled as a random 
process with memory. Unknown fading coefficients of the 
channels are estimated at the receivers with the assistance 



of the pilot symbols. We consider two channel estimation 
methods: single-pilot minimum mean-square-error (MMSE) 
estimation and noncausal Wiener filter estimation. We study 
AF and DF relaying techniques with two different transmission 
protocols. We obtain achievable rate expressions and optimize 
the training parameters by maximizing these expressions. 
We concentrate on the Gauss-Markov and lowpass fading 
processes for numerical analysis. 

II. Channel Model 

We consider a three-node-relay network which consists 
of one source node, one relay node and one destination 
node. Source-destination, source-relay and relay-destination 
channels are modeled as Rayleigh fading channels with fading 
coefficients denoted by hsd ^ CAf{0, fi^^), hsr ^ CAf{0, cr^^) 
and hrd ^ CAf{0, cr^jj^ respectively. Each channel is inde- 
pendent of others and exhibits memory with an arbitrary corre- 
lation structure. Hence {hsd}, {hsr}, and {hrd} are assumed 
to be mutually independent Gaussian random processes with 
power spectral densities Sh^^e^''^), Sh^^e^'^) and Sh^^ie'""), 
respectively. In this relay network, information is sent from the 
source to the destination with the aid of the relay. Transmission 
takes place in two phases: network training phase and data 
transmission phase. Over a duration of M symbols, the source 
and the relay are subject to the following power constraints: 

\\Xs4^+E{\\Xsf}<MPs (1) 
\\Xr4^ + +E{\\Xr\\^}<MPr (2) 

where Xg.t and x^.f are the source and relay pilot vectors, 
respectively, and and x^ are the data vectors sent by the 
source and the relay, respectively. 

III. Network Training Phase 

In the network training phase, source and relay send pilot 
symbols in nonoverlapping intervals with a period of M 
symbols to facilitate channel estimation at the receivers. In a 
block of Al symbols, transmission takes place in the following 
order. First, the source sends a single pilot symbol Xs,t, and 
the relay and destination receives 

yr,t = hsrXs,t + n-r and yd,t = hsdXs.t + rid, (3) 

~ CJ\f{m, (7^) is used to denote that x is a proper complex Gaussian 
random variable with mean m and variance . 



and estimates hgr and hsd, respectively. Then, transmission 
enters the data transmission phase, and source sends an 
{M — 2)/2-dimensional data vector that is again received 
by the relay and destination terminals. Next, only the relay 
sends a single pilot symbol Xr,t, and the signal received at the 
destination node is 

yd,t = hrdXr.t + n'd, (4) 

which is used by the destination to estimate hrd- In (O and (|4]l, 
fir, Ud, and are assumed to be independent and identically 
distributed (i.i.d.) zero mean Gaussian random variables with 
variance ct^, modeling the additive thermal noise present at the 
receivers. In the remaining duration of {M — 2)/2 symbols, 
transmission again enters the data transmission phase. In this 
case, the relay transmits an (A/ — 2)/2 dimensional data 
vector to the destination while the source either becomes silent 
or continues its transmission depending on the cooperation 
protocol. This order of transmission is repeated for the next 
block of M symbols. 

As noted before, we consider two channel estimation meth- 
ods. In the first method, only a single pilot symbol is used to 
obtain the MMSE estimate of the channel fading coefficients. 
As described in [5], MMSE estimates of the fading coefficients 
and the variances of the estimate errors are given as followS 
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where P^^ ^ and P^^ ^ are the power of the pilot symbols 
sent by the source and the relay, respectively, and y^.t ^ 
CAfiO,alP,^^^ +al), yd,t - CA/-(0, a^.F,^,, + and 

In the second method, we employ the noncausal Wiener 
filter which is the optimum linear estimator in the mean-square 
sense. The Wiener filter is employed at both the relay and 
the destination. Note that since pilot symbols are sent with a 
period of A/ symbols, the channels are sampled every MTg 
seconds, where Tg is the sampling time. As described in [10], 
we have to consider the undersampled versions of the Doppler 
spectrums of the fading coefficients, which are given by 
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and h are used to denote the estimate^and error in the estimate of h, 
respectively. Hence, we can write h = h + h. 



Then, the channel MMSE variances for the noncausal Wiener 
filter at time Ml + m are given by [11] 
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After obtaining the estimates, we can express the fading 
coefficients as 

hsr{Ml + to) = lisriMl + m) + hsr{Ml + to) 
hsd{Ml + to) = hsd{Ml + m) + hsd{Ml + to) 
hrd{Ml + m) = hrd{Ml + to) + hrd{Ml + m). (14) 

IV. Data Transmission Phase 

Note that as described in the previous section, within a block 
of M symbols, two symbol durations are allocated for channel 
training while data transmission is performed in the remaining 
portion of the time. We assume that relay operates in half- 
duplex mode. Hence, the relay first listens and then transmits 
to the destination. We consider two transmission protocols: 
non-overlapped and overlapped transmissions. 

A. Non-overlapped Case 

In this protocol, the source and relay send data symbols in 
nonoverlapping intervals. The source, after sending the pilot 
symbol, sends its (A/ — 2)/2 data symbols which are received 
by the relay and the destination a^ 



2/r,d(™) = hsr{m)xs4{m) + nr{m) 
Vd,d{f^) = hsd{m)xs,dim) + nd{m) m 
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Next, the source stops transmission, and the relay sends first 
its pilot symbol and then (A/ — 2)/2 data symbols which 
are generated from y^^ = [yr,d(2), . . . ,yrA{M /2)]. Thus the 
destination receives 

yd,d{j) = hrd{j)XrAj) + "dO') j = m + Af/2 (16) 

where again m = 2, . . . , Af/2. After substituting ( fT4] i into (fTSl l 
and ( fTSl l, we obtain 

Vr,d{'^) = hsr{m)xs4{m) + hsr{m)xs4{m) + nr{m) (17) 
VdA"^) = hsd{m)xs,d{m) + hsd{m)xs,d{m) + nd{m) 

yd,dU) = hrd{j)Xr,d{i) + hrd{j)Xr,d{j) + '^dU) 

where to = 2, ... , Ar/2 and j ^ m + M /2. 

^ Since we consider transmission in a block of M symbols, we drop the 
block index for the sake of simplicity and use m instead of using Ml + m. 
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B. Overlapped Case 

In this protocol, the source continues its transmission while 
the relay is sending its data symbols. The source becomes 
silent only when the relay is sending the pilot symbol. There- 
fore, the received signals in the data transmission phase can 
be written as 



yr,<i(m) = hsr{m)XsA{fn) + nr{m) 
yd,di'^) = hsdim)xsMim) +nd{m) 

yd,dU) = hsd{j)XsA3) + hrd{i)XrAi) + "-dO') 



(18) 



where m = 2, M/2 and j — m + M/2. Similarly as in the 
non-overlapped case, we can integrate the estimation results 
to (O and write 

yr,d(w) = hsrim)xsA{m) + hsrim)xsM{m) + ririm) (19) 
VdA^) = hsd{'m)xsA'^) + hsd{'m)xsAm) + nd{m) 
yd,dU) = hsd{j)xsA{j) + hsd{j)xsAj) 

+ Kd{j)Xr,dU) + Kd{j)Xr,d{j) + '^dU)- 

V. Achievable Rates for AF Scheme 

In this section, we consider the AF relaying scheme in 
which the relay sends to the destination simply the scaled 
version of the signal received from the source. An achievable 
rate expression for the AF scheme is obtained by maximizing 
the mutual information between the transmitted signal vector 
Xs d and the (M— 2)-dimensional received signal y^; ^ — 
[yd,d(2), . . . , ydAM/2), ydA^If^ + 2)^. • ■ ,ydAM)l given 
the estimates of the fading coefficients, h^^, i^sd and hrd are 
used to denote the vectors of channel estimates. Therefore, an 
achievable rate expression is given by 



Iaf = sup ^I{Xs^d;yd,d\^sr,'^sd,Kd)- 



(20) 



Note that the above formulation supposes that the destination 
node also knows hgr- Hence, it is assumed that these estimates 
are reliably forwarded by the relay to the destination using low 
rate links. A lower bound on Iaf can be obtained by assuming 
similarly as in [3] that the estimation errors are additional 
sources of worst-case Gaussian noise. We define the new noise 
random variables in non-overlapped and overlapped cases as 



ZrA"^) = hsr{m)XsA'^) + »^r(™) (21) 

ZdA^) = hsd{m)xsA^) + nd{m) 

ZdAj) = Kd{3)XrAj) + 
ZrA"^) = hsr{m)XsA'^) + "r(™) (22) 

ZdA'^) = hsd{m)xsA'>TT') + ridim) 

ZdAj) = hsd{i)XsAi) + hrd{j)XrAj) + "dO') 
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respectively. By assuming that the new noise components 
are Gaussian random variables and using techniques similar 



to those in [12], we can obtain the following worst-case 
achievable rate expression for the non-overlapped case: 
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and wsd ~ CM{0, 1), Wsr C7V(0, 1), Wrd CJ\f{Q, 1). 
Px^ dim) and Pj.^, ^(j) are the powers of the m*'* source symbol 
and j*'* relay symbol, respectively, and ct 
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Finally, note that j — m + M/2 

Similarly, we can find the following achievable rate expres 
sion for the overlapped case: 
1 

M 



(25) 
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and g(a, 6, c, d) = ^^'"^f^l'^"^ and j = m + A//2. 
VI. Achievable Rates for DF Scheme 

The repetition coding and the parallel coding are two 
possible coding techniques used in DF schemes [8]. First, we 
consider the repetition coding, and for this case the achievable 
rate is given by 

^rc = Y7 sup jmin /(x^^d; yr,d|hsr), /(xs,d; yd,d|hsd, hrd) } 

(27) 

Employing the techniques used in the AF non-overlapped 
scheme, we obtain the following achievable rate expression 
for non-overlapped DF with repetition coding: 



where 



nonov&r ,rc 



(28) 



h = log [1 + bi^A , h = log [1 + ai,„i + ci J-] , 

and ai.,„,6i^TO and ci.j are given in (l24l i. For the overlapped 
case of the DF repetition coding, (l28l l holds with Ii and I2 
defined as 

h = log [1 + C2j] , h = log [1 + a2,m + 62J + C?2,m + a2,m&2j] 
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Fig. 1. The optimal achievable rates vs. SNR for the Gauss-Markov fading 
model (a = 0.99) and different relaying techniques, a^^ = l,crj^ = 16 and 
(T^^ = 16. (S: single-pilot estimation. W: Wiener filter.) 



Fig. 3. The optimal achievable rates vs. SNR for the lowpass fading model 
when noncausal Wiener filter is employed, a^^ = 1, cr^^ = 4 and ct^^ = 4. 
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Fig. 2. The optimal achievable rates vs. SNR for the lowpass fading model 
when noncausal Wiener filter is employed, a^^ = 1, cr^^ = 16 and cr^^ = 16. 



Fig. 4. Normalized bit energies -j^ vs. SNR for the lowpass fading model 



when noncausal Wiener filter is employed, a , = 1, tr^ 



16 and a^, 

rd 
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where a2,m,b2.j,C2.j, and (i2,m are given in (1261 1. When we 
employ the parallel coding, we have 

Ipc = Y7 ^29) 

|min /(x^,d;y^,d|hsr),/(xs,d;yd,d|h.d) +^(xr,d;y<j,d|hrd) }• 

Similarly, we can find, for the nonoverlapping case, an achiev- 
able rate expression given by 



= -^E.^s.En,^^E^^^ ^ min(/i, /2) 



where 

Ji =log(l + 6i,„0, /2 = log(l 



(30) 



l0g(l + Ci,j). 



VII. Optimizing Training Parameters 

In this section, we consider two particular fading processes. 
In the first case, fading is modeled as a first-order Gauss- 
Markov process whose dynamics is described by 

hk — ahk-i + Zk < a < 1 fc = 1, 2, 3, . . . 

where {zk} are i.i.d. circular complex Gaussian variables with 
zero mean and variance equal to (1-Q!^)(t^. In the above 



formulation, a is a parameter that controls the rate of the 
channel variations between consecutive transmissions. The 
power spectral density of the Gauss-Markov process with 
variance cr^ is given by 



(31) 



1 + — 2acos{w) 

We also model the fading as a lowpass Gaussian process whose 
power spectral density is given by 

S.ien - ( W < (32) 

[ 0, otherwise 

where Wd — "^Trfd is the maximum Doppler spread in radians. 

In Gauss-Markov channels, it is difficult to find a closed- 
form expression for the variance of the estimate error when 
Wiener filter is used, because the channel's spectrum is not 
band limited. Therefore, there is always aUasing in the un- 
dersampled Doppler spectrums, which causes an increase in 
the variance of the error. On the other hand, when fading is 
modeled as a lowpass process, we can find a explicit solution 
for the error variance, and we can express it as 
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Fig. 5. Normalized bit energies vs. SNR for the lowpass fading model 
when noncausal Wiener filter is employed, u 



■ 4 and cr^ , : 
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In the lowpass case, if the channel is sampled sufficiently 
fast (i.e., M < tjV), there is no aliasing and the power is 
distributed equally among data symbols. However, note that 
the power allocated to the data symbols of the source is not 
equal to the power allocated to the data symbols of the relay. 
In general, if there is aliasing or a single pilot is used for 
estimation, the power allocated to the data symbols will differ 
depending on their distance to the pilot signals. 

Having obtained achievable rate expressions, our next goal 
is to jointly optimize training period M, training power, and 
power allocated to the data symbols. 

VIII. Numerical Results 

In this section, we present numerical optimization results. 
In Figure [T] we plot the optimal achievable rates with respect 
to SNR for different relaying protocols by using two different 
methods of channel estimation. Fading is assumed to be a 
Gauss-Markov process. The dashed lines indicate the optimal 
data rates obtained when noncausal Wiener filter is used, 
whereas the solid lines show the optimal data rates obtained 
when a single pilot symbol is used for estimation. The rates 
are optimal in the sense that they are obtained with optimal 
training parameters and optimal power allocations. We can 
see that at low SNR values, DF provides higher rates and 
parallel non-overlapped DF scheme is the most efficient one. 
As expected, Wiener filter performance is better than that 
of the estimation that uses a single pilot. Moreover, at low 
SNR values non-overlapped and overlapped relaying schemes 
give the same optimal results, and optimal power distributions 
among data and pilot symbols are the same for both. On the 
other hand, at high SNR values, we see a significant increase 
in the data rate of AF overlapped scheme compared to the 
other schemes. 

In Fig. |2] and Fig. |3] we plot the optimal data rates when we 
estimate the lowpass fading process using a noncausal Wiener 
filter. The channel variances are cr^^ — 1, CTj^ — 16, a^^ = 16, 
and cr^^ = 1, CTg^ = 4 and a^^ — 4, respectively. Conclusions 
similar to that given for Fig. [T] are drawn again. In Figs. |4] 
and|5] the bit energy normalized by the noise variance, is 
plotted as a function of SNR. In all cases, we observe that 



Fig. 6. The optimal acliievable rates vs. training period M for tlie GausS' 
Markov fading model. Single-pilot MMSE estimation is employed. The 
dashed lines are obtained when a^, = 1, (Tg, 
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7. Optimal power distribution among the pilot and data symbols when 
5^ = 16, cr^^ = 16, SNR=OdB. Fading is a Gauss-Markov process 
0.99. Single-pilot MMSE estimation is employed. The optimal 
period is M = 16. Note that the first 8 symbols belong to the source and the 
last 8 bars belong to the relay. 
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Fig. 8. Optimal power distribution among the pilot and data symbols when 
(T^^ = l,o-g^ = 4,(T^^ = 4, SNR=OdB. Fading is a Gauss-Markov process 
with a = 0.99. Single-pilot MMSE estimation is employed. The optimal 
period is M = 30. Note that the first 15 symbols belong to the source and 
the last 15 bars belong to the relay 
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Fig. 9. Optimal power distribution among the pilot and data symbols when 
(j^^ = IjC^^ = 16, (J^^ = 16, SNR=OdB. Fading is a Gauss-Markov process 
with a = 0.99. Wiener filter is employed. The optimal period is M = 12. 
Note that the first 6 symbols belong to the source and the last 6 bars belong 
to the relay 
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Fig. 10. Optimal power distribution among the pilot and data symbols when 
(T^^ = IjCj^ = 16, (Tj;^ = 16,SNR=0dB. Fading is a Gauss-markov process 
with a = 0.99. Single-pilot MMSE estimation is employed. The optimal 
period AI = 12. Note that the first 6 symbols belong to the source and the 
last 6 bars belong to the relay 



minimum bit energy is achieved at a nonzero SNR value. 
If SNR is further decreased, higher bit energy values are 
required and hence, operation at these very low SNRs should 
be avoided. 

In Figure |6] we plot the optimal data rates as a function of 
the training period, M, when SNR=0 dB for different relaying 
schemes and different channel variances. Single-pilot-symbol 
estimation is employed. Since a relatively low SNR value is 
considered, AF non-overlapped and AF overlapped schemes 
provide lowest rates. The highest performance is obtained 
when DF parallel non-overlapped scheme is used. 

In Figs. |7] and [8] power allocated to the pilot and data 
symbols is plotted when Gauss-Markov channel is considered 
and AF non-overlapped scheme is employed. The first half 
of the bars shows the power allocated to the source symbols 
and the rest shows the power allocated to the relay symbols. 
The first bar of the each group gives the power of the pilot 
symbols. Note that these power distributions are obtained 
when the period is at its optimal value when SNR=0 dB. 



The optimal periods are 16 and 30 when a^^ — l,cr^^ — 16, 
^rd — 16, and cr^^ ~ ^-.(^Ir — 4, ~ 4, respectively. In 
Figure |9] the optimal power distribution is displayed when 
noncausal Wiener filter is used for estimation at SNR = 
dB. Note that the optimal period is 12. In Figure [TOl we 
plot the power distribution when single-pilot estimation is 
performed at the optimal period 12 of Wiener Filter estimation. 
It is observed that more power is given to the pilot symbols 
when single-pilot-symbol estimation is used. However, when 
we employ noncausal Wiener filter, the power allocated to the 
pilot symbols is decreased thereby increasing the data rate by 
giving more power to the data symbols. 

IX. Conclusion 

We have studied transmission over imperfectly-known relay 
channels. The channels are learned using single-pilot MMSE 
estimation or noncausal Wiener filter. We have obtained 
achievable rate expressions for both AF and DF relaying 
schemes. Subsequently, we have jointly optimized the training 
period and power, and data power levels in Gauss-Markov and 
lowpass fading models. We have compared the performances 
of different relaying techniques at different SNR values and 
different channel variances. 
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